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Abstract. The behaviour, ecology, and phylogenetics of glassfrogs from the Atlantic Forest are poorly understood. To con- 
tribute to these three topics, we studied reproductive behaviours, microhabitat use, influence of environmental variables 
on the activity of adults, and clutch thickness, as well as the phylogenetic relationships of V. baliomma, a key centrolenid 
species from the northern Atlantic Forest. We sampled two streams in the south of the state of Bahia, Brazil, between No- 
vember 2018 and April 2019, and measured relevant specimen and environmental variables. Our results include (i) the first 
observations of female-only egg-brooding in this species, (ii) the statistically significant dependence of reproductive activ- 
ity on daily rainfall and high humidity, and (iii) that adults use different leaf face for clutch deposition and activity, with 
clutches being placed on both sides of leaves and at statistically significant lower heights and horizontal distances from 
water. Maximum likelihood phylogenetic analyses of DNA sequences of three mitochondrial and one nuclear gene (12, 16S, 
COI, and POMC, respectively) confirmed the phylogenetic position of V. baliomma within the genus. Although our results 
inferred V. baliomma as sister of V. gorzulae (a Guiana Shield species) and both as sister taxa of all other Atlantic Forest and 
Cerrado species, these relationships have low support. 
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Introduction 


Anurans have a large diversity of traits that are related to 
their reproductive biology (DUELLMAN & TRUEB 1994) 
and are one of the groups with the greatest variety of re- 
productive modes among vertebrates (HADDAD & PRADO 
2005, CRUMP 2015). Their reproductive aspects provide im- 
portant insights into their systematics, evolution, ecology, 
and may determine species conservation measures (HOFF- 
MANN et al. 2008, GUAYASAMIN et al. 2009, VARGAS-SALI- 
NAS et al. 2014, DELIA et al. 2017, ESCALONA-SULBARAN et 
al. 2019). However, information on reproductive biology is 
scarce or even nonexistent for many species (WELLS 2007, 
DELIA et al. 2017). 

The family Centrolenidae is an interesting group due to 
its morphological, ecological, and behavioural character- 
istics (RUIZ-CARRANZA & LYNCH 1991, SENARIS & AYAR- 
ZAGUENA 2005, CISNEROS-HEREDIA & MCDIARMID 2007), 
as well as its phylogenetic and biogeographical complex- 
ity (GUAYASAMIN et al. 2008, CASTROVIEJO-FISHER et 


al. 2014). This family is restricted to the American Trop- 
ics, where its species inhabit forest-lined waterfalls and 
streams, perching on vegetation and laying their eggs on 
leaves, branches, mosses, or rocks overhanging the water 
(CISNEROS-HEREDIA & MCDIARMID 2007, GUAYASAMIN 
et al. 2009, Twomey et al. 2014). After hatching, tadpoles 
drop into the water and continue their development buried 
in sand and debris on the bottom of water bodies (KUBICKI 
2007, GUAYASAMIN et al. 2009). 

While most of the few described species were known 
from cursory descriptions of their external morphology 
based on preserved material before the 1990’s (e.g., TAY- 
LOR 1949, RIVERO 1968, LYNCH & DUELLMAN 1973, DUELL- 
MAN 1976, RIVERO 1985), centrolenids are now emerging as 
a model system for ecological and evolutionary studies (e.g., 
HUTTER et al. 2013, DELIA et al. 2013, 2017, 2020, CASTRO- 
VIEJO-FISHER et al. 2014, ESCALONA-SULBARAN et al. 2019, 
BARNETT ET AL 2020). In fact, glassfrogs have become a re- 
current model for studies addressing amphibian reproduc- 
tive strategies, with focus on issues such as parental care 
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(e.g., HAYES 1991, VOCKENHUBER et al. 2009, DELIA et al. 
2017, 2020, SCHULTE et al. 2020), acoustic communication 
(e.g., WELLS & SCHWARTZ 1982, ESCALONA-SULBARAN et 
al. 2019), reproductive activity and patterns (e.g., GREER & 
WELLS 1980, JACOBSON 1985, BasTO-Riascos et al. 2017, VA- 
LENCIA-AGUILAR et al. 2020a), and microhabitat selection 
(e.g., VALENCIA-AGUILAR et al. 2012, CABANZO-OLARTE et 
al. 2013, DÍAZ-RICAURTE et al. 2019, OrTIZz-Ross et al. 2020). 

Parental care in glassfrogs was considered a rare be- 
haviour, exclusive to males (MCDIARMID 1978), but recent 
studies reported that this is in fact a widespread trait in 
this family, with an interesting evolutionary pattern in var- 
iation of time invested (brief or prolonged) and sex (male 
or female) of providing parental care (BRAVO- VALENCIA 
& DELIA 2016, DELIA et al. 2017, 2020). Parental care in- 
creases the chances of survival or benefits the growth of 
offspring, thus being an important measure to maximize 
fitness (TRIVERS 1972, MCDIARMID 1978, CLUTTON-BROCK 
1991). Brooding behaviour is one mode of egg care found 
in the Centrolenidae. It is performed by ventral contact, 
where the parent positions its body over the clutch (HAYES 
1991, GUAYASAMIN et al. 2009) in a quest to protect em- 
bryos from desiccation, predation, developmental abnor- 
malities, and fungal infection (HAYES 1991, CRUMP 1995, 
VOCKENHUBER et al. 2009, DELIA et al. 2013, 2017, LEHTI- 
NEN et al. 2014, BRAVO- VALENCIA & DELIA 2016). 

Reproductive site selection (call and oviposition site) is 
another important trait related to parental fitness and off- 
spring survival (RESETARITS 1996). Reproductive micro- 
habitat has also been suggested to influence centrolenid 
call variability (BoscH & DE LA Riva 2004) - for exam- 
ple, it increases the effectiveness of acoustic signalling in 
Hyalinobatrachium fleischmanni (GREER & WELLS 1980, 
WELLS & SCHWARTZ 1982) - and apparently plays an im- 
portant role in the evolution of anuran sexual size dimor- 
phism (Sitva et al. 2020). Furthermore, the oviposition 
site is of great importance to offspring survivorship (RESE- 
TARITS 1996). Females make choices from a complex com- 
position of options that vary according to both biotic and 
abiotic cues, such as absence of predators and desiccation 
risk (BERNARDO 1996, RUDOLF & RODEL 2005, TOUCHON 
& WORLEY 2015). Oviposition site choice in centrolenids 
appears to protect embryos from the very same risks that 
favour parental attendance of eggs, with an evolutionary 
association between time spent on care, clutch structure, 
and oviposition site (DELIA et al. 2020). 

Within glassfrogs, the genus Vitreorana, with ten spe- 
cies distributed throughout the Cordillera de la Costa, 
Guiana Shield, Amazonia, Atlantic Forest, and Cerrado, 
has received less attention than other genera (GUAYASAMIN 
et al. 2009, PONTES et al. 2014, SANTANA et al. 2015, FROST 
2020). For example, Vitreorana baliomma is a small species 
(snout-vent length = 16.4-21.2 mm) described only based 
on morphological characteristics of preserved specimens 
(PonTES et al. 2014), and with its advertisement call quan- 
titatively described just recently (BANG et al. 2020). It oc- 
curs in forests associated with streams and creeks in the 
northern portion of the Atlantic Forest of Brazil, from 100 
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to 600 m a.s.l. (GOUVEIA et al. 2012, PONTES et al. 2014). 
Individuals use marginal vegetation along streams, perch- 
ing on leaves and trunks, ranging in height from 0.3 to 4m 
(GouVEIA et al. 2012). GOUVEIA et al. (2012), and LISBOA 
et al. (2019) presented notes about the reproduction of 
V. baliomma, indicating that females lay their eggs, pref- 
erably, but not exclusively, on the upper side of leaves, at 
0.5 to 2.0 m high, although neither study detected signs of 
parental care. 

Other fundamental aspects of its reproduction, such 
as microhabitat use or whether climatic conditions influ- 
ence behavioural activity, are also unknown. This lack of 
data is problematic if we consider the phylogenetic po- 
sition of Vitreorana as one of the earliest splitting clades 
within the tribe Cochranellini (GUAYASAMIN et al. 2008, 
CASTROVIEJO-FISHER et al. 2014, Twomey et al. 2014) and 
that the aforementioned characters are also unknown for 
most of the other Vitreorana species. For example, V. ri- 
tae is one of the species for which maternal care has been 
reported (formerly no care), nevertheless, just one speci- 
men was observed (DELIA et al. 2017), and we do not as yet 
know whether this behaviour is possibly more widespread 
in this genus. 

The phylogenetic affinities of V. baliomma are specu- 
lative and even its placement within Vitreorana is some- 
what tentative, being based upon the following morpho- 
logical characters: parietal and urinary bladder peritonea 
transparent, pericardium and visceral peritoneum white, 
bones green in life, liver bulbous, and dorsal coloration in 
preservative lavender (PonTEs et al. 2014). Besides Vitreo- 
rana, this combination of characters is also present, with- 
in the Centrolenidae, in species of the genera Chimerella, 
Hyalinobatrachium, and Teratohyla. Faced with this situa- 
tion and in the absence of DNA sequence analyses, Pon- 
TES et al. (2014) opted for allocating this species to Vitreo- 
rana. Their decision was based on biogeographic affinities, 
because all known species of glassfrogs from the Atlantic 
Forest form a clade within Vitreorana. This Atlantic For- 
est clade is the sister taxon of the Amazonian and Guiana 
Shield species of the genus (GUAYASAMIN et al. 2008, 2020, 
CASTROVIEJO-FISHER et al. 2014, Twomey et al. 2014). 
However, all the Atlantic Forest species of Vitreorana have 
the synapomorphy of a white urinary bladder peritoneum 
(GUAYASAMIN et al. 2009), which is conspicuously absent 
in all other species of Vitreorana (including V. baliomma), 
and casts doubts on the phylogenetic affinities of V. bali- 
omma. The phylogenetic relationships of V. baliomma 
transcend the details of glassfrog systematics because the 
relationships between the Atlantic Forest and Amazonia- 
Guiana Shield clades constitute a key example of the bio- 
geographic connection between these two regions in the 
late Miocene (e.g., CASTROVIEJO-FISHER et al. 2014, ROD- 
RIGUES et al. 2014, PRATES et al. 2017). 

In this study, we identified Vitreorana baliomma as a 
species of primary interest to overcome important knowl- 
edge gaps in the phylogenetic relationships (Darwinian 
shortfall) and biological traits (Raunkizran shortfall) of 
glassfrogs (HorTAL et al. 2015). Although modern tech- 
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nology may offer partial solutions to these shortfalls, field- 
work is still the most important means to fill these gaps 
(WILSON 2017). Considering the aforementioned relevance 
of V. baliomma, we conducted field expeditions to collect 
data to address the following issues: (1) describe the pattern 
of reproductive behaviours; (2) investigate the possible oc- 
currence of parental care; (3) characterize the physical pa- 
rameters of the sites used by adults for activity and deposit- 
ing egg clutches; (4) characterize egg clutches; (5) evaluate 
the influence of abiotic variables on the number of active 
adult individuals; (6) evaluate the influence of microhabi- 
tat on clutch thickness; and (7) infer its phylogenetic rela- 
tionships within Centrolenidae using DNA sequences and 
quantitative phylogenetic analyses. 


Material and methods 
Study sites 


We collected data between November 2018 and April 2019, 
in two slow- to moderate-current streams with sandy/ 
rocky beds: (1) Fazenda Santo Antônio (15°30°8.65” S, 
39°22'33.14” W; 105 m a.s.l.), municipality of Mascote and (2) 
at 3 km of Cachoeira do Tijuipe (14°23 17.61” S, 39°2’37.89” W; 
89 m a.s.l.) near BAoot1 - Ilhéus-Itacaré road, municipality of 
Urucuca (hereafter Mascote and BAoo1, respectively). Both 
streams are located in Atlantic Forest in the south of the 
state of Bahia, Brazil. The stream in Mascote runs through a 
cabruca forest, which is an agroforestry system where some 
of the native forest serves to provide shade for cocoa trees 
(Theobroma cocoa). The BAoo1 locality constitutes a well- 
regenerated secondary forest. The sampled area lies in the 
conservation unit Area de Proteção Ambiental (APA) Costa 
de Itacaré and within the buffer zone of the Parque Estadual 
Serra do Conduru (PESC). We collected all the data used for 
our reproductive behaviour studies from a 150-m transect in 
Mascote and a 100-m one in BAoo1, where we observed the 
greatest aggregation of individuals. In southern Bahia, the 
climate is characterized as being hot and humid, without a 
dry season (KOPPEN 1936). Mean annual rainfall is between 
1500-1700 mm, mean annual temperature is 23-24°C, and 
relative humidity is > 80% (ROEDER 1975). To better un- 
derstand the geographic distribution of the species, we re- 
viewed the literature and the Museu de Zoologia da Univer- 
sidade Estadual de Santa Cruz (MZUESC) scientific collec- 
tion from which we here provide an updated distribution 
map for this species (Fig. 1). 


Fieldwork 


We sampled at night, between 17:00 and 03:00 h, with spo- 
radic diurnal visits to verify the permanence of adults near 
clutches, which totalled 32 hours of sampling within 13 
days. We measured temperature and humidity at the begin- 
ning and end of each sampling period (period of time from 
arrival until leaving the stream), using a digital thermo- 
hygrometer Instrutherm HT-260 (accuracy + 2.5°C; + 3.5% 


RH). For each sampling period, we used as a descriptive 
variable the sum of the start and end values divided by 
two. Also, we measured daily precipitation for Mascote 
using a 7o-mm pluviometer fixed at 5 cm atop a wooden 
stake, placed in a locale without canopy cover. We installed 
the pluviometer just before sampling (between 16:30 and 
16:50 h) and read it every 24 hours. For BAoo1, we used 
daily precipitation values (00:00-01:00 to 23:00-00:00 h) 
from the Instituto Nacional de Meteorologia (INMET - 
https://portal.inmet.gov.br/), as recorded by the Ilhéus au- 
tomated weather station (14°39731.69” S, 39°10°53.18” W; 
80 m a.s.l.), located 33 km from our sampling site. 

We conducted visual encounter surveys for adults 
and clutches and auditory searches for calling males. We 
marked encounter sites of adults and clutches with unique- 
ly coded flagging tape. When possible, we differentiated 
males from females by observing the vocal sac, present 
only in males, and by the presence of oocytes in the ovi- 
duct, visible through the translucent belly of reproductive- 
ly active females. We collected at random representative 
voucher specimens of adults (females and males), embry- 
os, and tadpoles. We euthanized specimens in the lab with 
a topical dose of the anaesthetic benzocaine (50 mg/g) and 
fixed them in 10% formalin. After 24 h, we preserved adults 
in 70% ethanol. We removed samples of thigh muscle tis- 
sue from some specimens and preserved them in 96% eth- 
anol for DNA analyses. We deposited specimens and tis- 
sue samples in the Amphibian Collection of the Museu de 
Ciências e Tecnologia da Pontificia Universidade Católica 
do Rio Grande do Sul (MCP). 


Behavioural observations 


We opportunistically monitored pairs in amplexus until 
oviposition and recorded other intraspecific interaction 
events throughout the duration of amplexus. We used a 
headlamp with a red-light filter to minimize disturbance 
(VaRGAS-SALINAS et al. 2014, ToUZOT et al. 2020). We vis- 
ited egg clutches during the following day and night to 
evaluate the presence of adults and their potential interac- 
tions with egg clutches after oviposition. We recorded the 
number of males heard vocalizing and the number of indi- 
viduals visually detected (males and females) during each 
sampling period to assess if the measured abiotic variables 
influenced the species’ activity (see Statistical analyses). We 
recorded these data across the 13 days of sampling, with a 
frequency of at least one-night monthly in one of the lo- 
calities and duration of sampling periods between 2-5 h. 


Microhabitat data 


We recorded the following variables for each encounter 
site of adults and clutches: (1) vertical distance to substrate, 
(2) horizontal distance from streambed, (3) type of struc- 
ture used as perch (i.e., leaf, branch, moss, rock), (4) perch 
size, (5) face (i.e., upper or lower), and (6) type of leaf sur- 
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face (i.e., smooth or hairy), (7) position relative to the stre- 
ambed (i.e., streambed, bank of stream, outside the stream 
flood area), (8) side to the stream flow (i.e., left or right), 
and (9) canopy cover. We measured the height of the ex- 
act encounter site above the substrate, and the horizontal 
distance of the location of the adult, pair, or clutch on the 
perch to the streambed margin. We quantified both varia- 
bles using a 5-m millimetric tape. We quantified perch size 
using a 2-m millimetric tape to measure the length (L) and 
width (W) of the leaf and calculated the elliptic area (A) 
with the formula 


A=kx Wx on, 


We quantified canopy cover using photos taken during 
daytime visits, with a fisheye lens attached to a smartphone 
camera. We took all photos with the phone positioned 5 cm 
above the centre of the encounter site and levelled at 180°. 
For each photo, we calculated the percentage of canopy- 
free pixels using the software Gap Light Analyzer (FRa- 


ZER et al. 1999). We compared adult and clutch data to test 
if there is a difference between the microhabitats used by 
adults and oviposition sites (see Statistical analyses). 


Egg clutch characterization 


We recorded the diameter and thickness of the entire egg 
capsule, the number of eggs, stage of development, and 
embryo coloration of each observed clutch. We measured 
length (L) and width (W) using a digital calliper (accuracy 
+ 0.02 mm), and then calculated the diameter (D) by us- 
ing the formula 


= A 
D=2 T, 


by which we defined elliptic area (A) for clutches, too. To 
measure clutch thickness, we used small wooden rods 
scaled in millimetres (mm). We used GOSNER (1960) to 
identify embryonic stages. Moreover, we tested if the char- 
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Figure 1. Map of northeastern Brazil with symbols indicating known localities of Vitreorana baliomma. 1 - Fazenda Novo Pau Brasil, 
municipality of Itamaraju; 2 - Fazenda Santo Antônio, municipality of Mascote; 3 - BA001 (Ilhéus-Itacaré road), municipality of 
Uruçuca; 4 - RPPN Ararauna, municipality of Una; 5 - Parque Estadual Serra do Conduru, municipality of Uruçuca; 6 - Parque 
Nacional (PARNA) da Serra das Lontras, municipality of Arataca; 7 - PARNA Serra de Itabaiana, municipality of Areia Branca; 
8 - Fazenda Palmeiras, municipality of Itapebi; 9 - Fazenda Limoeiro, municipality of Almenara; 10 - Murici Ecological Station, 
municipality of Murici. Brazilian states included are: Alagoas (AL), Bahia (BA), Minas Gerais (MG) and Sergipe (SE). We collected 
reproductive behaviour data from localities 2 and 3, while sequenced tissues are from 2, 3, 4 and 5. 
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acteristics of clutch structure or oviposition site are related 
to clutch thickness (see Statistical analyses). We assigned 
clutches to V. baliomma based on the following criteria and 
observations: 1) clutches, egg capsules, and embryos were 
similar in structure, colour, size, and position site to those 
laid by females directly observed by us; 2) despite our field 
efforts, we did not find or hear other species of glassfrogs 
in the area around the sampled streams; 3) egg clutches of 
other known sympatric anuran species that lay their eggs 
on leaves (e.g., Dendropsophus spp. and Phyllomedusa spp.) 
are clearly different from those of V. baliomma. For ex- 
ample, females of Phyllomedusa spp. build a nest with the 
leaves; 4) we incubated one clutch in the laboratory and 
observed its development from oviposition to hatching, to 
identify the developmental stages. 


Molecular methods and phylogenetic analyses 


We extracted genomic DNA from ethanol-preserved 
muscle tissue using the Wizard Genomic DNA Purifica- 
tion Kit (Promega, Madison-WI, USA) as per manufac- 
turer’s instructions. We performed polymerase chain re- 
actions (PCR) in 25-uL volumes with 2.0 UL of genomic 
DNA to amplify ~ 646 bp of the COI mitochondrial gene 
(T3-AnF1 and T7-AnRi; Lyra et al. 2017), ~576 bp of the 
16S rRNA gene (primers 16Sar and 16Sbr; PALUMBI et al. 
1991), 562 bp of the 12S rRNA gene (primers t-Phe-frog 
and 12S-frogR; WIENS et al. 2005), and ~492 bp of POMC 
(primers POMC-1 and POMC-2; Wrens et al. 2005). We 
checked amplifications in agarose gels (1-2%), sequenced 
fragments in both directions, and used Sequencher v4.1.4 
(Gene Codes Co., Ann Arbor-MI, USA) to assemble and 
edit the resultant chromatograms. 

Additionally, we downloaded homologous sequenc- 
es of other Vitreorana spp. available from GenBank (as of 
30 December 2019). As outgroups, we used homologous se- 
quences from GenBank of one or two terminals of all other 
centrolenid genera. We rooted all trees with sequences of 
Allophryne ruthveni, which forms part of the sister group of 
Centrolenidae (GUAYASAMIN et al. 2008, 2009). For some 
terminals, if genetic distances were <1%, we assigned se- 
quences from different specimens to a single composite ter- 
minal (Supplementary Table S1). We did this to reduce the 
numbers of missing entries per terminal and search space 
during phylogenetic analyses (WILKINSON 1995, KEARNEY 
2002, BROWER 2018). Our final dataset included a total of 
21 species and 185 terminals (12 for outgroup and 173 for in- 
group). Specimens and DNA sequences used in the phylo- 
genetic analysis are listed in Supplementary Table S1. 

We individually aligned sequences of each marker using 
MUSCLE (EDGAR 2004) in AliView (LARSSON 2014) un- 
der default parameters. We combined the alignments men- 
tioned above into a single matrix using SequenceMatrix 1.8 
(Varpya et al. 2011). Phylogenetic analysis was based on 
the maximum likelihood (ML) criterion. We determined 
the most appropriate model of nucleotide substitution and 
data partition using ModelFinder (KALYAANAMOORTHY et 


al. 2017) in IQ-TREE. We divided the dataset by gene and, 
for COI and POMC, by codon positions. We performed 
500 independent searches in IQ-TREE (Minu et al. 2020) 
under default settings and using the models and partitions 
derived from the results of ModelFinder (CHERNOMOR et 
al. 2016). We assessed clade support with 1000 bootstrap 
(BS) pseudoreplicates using IQ-TREE with the same afore- 
mentioned settings. We used the SumTrees package in 
DendroPy v4.4.0 (SUKUMARAN & HOLDER 2010) to add BS 
frequencies to the best ML tree. 


Statistical analyses 


We performed all analyses in R 4.0.2 (R Core Team, 2020) 
and RStudio 1.3.1073 (RStudio Team, 2020) code edi- 
tor, with packages Ime4 (Bates et al. 2015) for the Linear 
Mixed Model (LMM) and Generalized Linear Mixed Mod- 
el (GLMM) analyses, car (Fox & WEISBERG 2019) to obtain 
the significance of models using likelihood ratio test, AIC- 
cmodavg (MAZEROLLE 2020) for model selection based on 
AICc and ggplot2 (WICKHAM 2016), and gridExtra (Au- 
GUIE 2017) to generate the figures. Other functions used 
are those present in the R base package. We considered re- 
sults as statistically significant if P-values were < 0.05. 

Prior to analysis, we scaled and centred (function: scale) 
temperature, humidity and precipitation to standardize ex- 
planatory variables and logarithm-transformed (function: 
log) leaf size values to decrease data amplitude. We ran 
Shapiro-Wilk tests (function: shapiro.test) to test the nor- 
mality of dependent variables for models and quantitative 
variables for comparison of microhabitat use. We found 
that the nightly number of calling males (W = 0.83706, 
p = 0.01945), sampled individuals (W = 0.66706, p = 
0.0002649), horizontal distance (W = 0.55034, p < 0.001), 
leaf size (W = 0.94141, p = 0.02917), and canopy openness 
(W = 0.73687, p < 0.001) all followed a non-normal distri- 
bution. Clutch thickness (W = 0.97026, p = 0.8026) and 
vertical height (W = 0.97734, p = 0.5457) instead followed 
a normal distribution. 

Consequently, we performed a GLMM analysis (func- 
tion: glmer) to test if average temperature and humidity, 
and daily rainfall are predictors of the number of calling 
males and sampled individuals. We used the sampling lo- 
cality as a random variable, while considering different 
abundances in each area, in both analyses. 

To test if the variables summarizing clutch structure 
(ie., number of eggs, diameter, and developmental stage) 
are predictors of clutch thickness, we performed a Line- 
ar Model (LM) analysis (function: Im). After that, we per- 
formed an LMM analysis (function: lmer) to test if micro- 
habitat variables are predictors of clutch thickness, and se- 
lected just vertical height, leaf size, and leaf surface as in- 
dependent variables to prevent overparameterized models 
based on our sample size. We added developmental stage 
as a random effect to test if microhabitat variables are pre- 
dictors of clutch thickness in LMM, given that it is an in- 
trinsic feature of clutches (see Results). 
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We selected additive or simple models using the Akaike 
information criterion (AIC; AKAIKE 1973). To generate 
AICc comparable values for models, we performed the 
function aictab. For each dependent variable, we fitted a 
simple model to each independent variable, an additive 
model with independent variables without interaction, 
and an additive model with independent variables inter- 
acting between them. After selection, we used the analysis 
of variance (ANOVA; function: Anova) as a likelihood ra- 
tio test to quantify the significance of variables in models. 
We fitted LMM and GLMM with the likelihood function 
and Poisson distribution, respectively. 

To assess if there were differences between the micro- 
habitat used by adults and oviposition sites, we performed 
independent 2-group Students t-Tests (function: t.test) for 
variables with normal distributions, and Mann-Whitney U 
Tests (function: wilcox.test) for variables with non-normal 
distributions, and beyond these, Pearson’s Chi-squared 
Tests (function: chisq.test) for categorical variables (i.e., 
leaf surface, position at stream, stream side). 


Results 


We found a total of 24 individuals (14 males, six females, 
and four unsexed individuals ) and 19 clutches (we disre- 
garded one clutch with hatched embryos for its number of 
eggs, diameter, and thickness). During sampling, tempera- 
ture varied from 23.7-29.1°C, humidity from 77-93.2%, and 
daily rainfall was between o and 53 mm. The region of the 
study area experienced abnormally low rainfall between 
January and March 2019 compared to the previous 10 years. 


Reproductive behaviour 


We observed four mating events, one in its entirety (in- 
cluding courtship, amplexus, oviposition, associated be- 
haviours, and female egg brooding), two without preced- 
ing courtship, and one without courtship, oviposition, and 
brooding. On the days of our observations, temperatures 
varied from 24.3-26°C, humidity from 90.7-93.2%, and 
daily rainfall from 14-18.4 mm. Our description is based 
mainly on a pair encountered on 11 April 2019, in which 
we could observe all behaviours from the beginning to the 
end. At approximately 22:15 h, we discovered a female po- 
sitioned ~ 30 cm away from a male vocalizing from a bush. 
The male was calling from the upper side of a leaf ~ 100 cm 
above the ground. The female moved across the foliage, cir- 
cling the male's calling site. At 22:25 h, the female leaped to 
the leaf positioned in front of the calling male. The male, 
still calling, turned to face the female, which then leaped 
on to the same leaf. The male then climbed on the female's 
back and initiated an axillary amplexus. The amplectant 
pair moved over the bush by the female leaping to leaves 
above and below the initial encounter site (Fig. 2A) and in 
the process switching between the upper and lower faces 
of the leaves (Fig. 2B). While on a leaf, the pair constant- 
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ly performed 180 and 360° rotations (Fig. 2C). After 102 
minutes in amplexus, the female deposited a clutch on the 
upper surface of a leaf. Right after oviposition, the female 
moved backwards, positioning her body above the clutch in 
a brooding posture. The male released the female just after 
she had laid her clutch, leaped to a leaf approximately 10 cm 
above the oviposition site, and began to vocalize again. The 
female remained positioned above the clutch for another 
110 minutes, raising the posterior part of her body at inter- 
vals of 1 to 10 minutes (Fig. 2D). At 01:58 h, the female left 
the clutch, remained next to it for a few seconds more, then 
leaped to other leaves and left the oviposition site. 

In the other three mating events, we observed the same 
behaviours as exhibited by the pair described above, with 
the following differences. In one case, at 21:10 h, also on 
11 April 2019, we found the individuals already in amplexus 
(continued for a further 25 min), with the female deposit- 
ing the clutch on the lower face of the leaf and brooding 
the eggs. In the third mating event, we found individuals 
in amplexus at 23:14 h, on 14 April 2019. While in amplex- 
us, the female performed movements with its front limbs, 
rubbing them on her face. After 56 min in amplexus, it 
started to rain and the pair became more active, moving 
about more frequently in the vegetation. At 01:08 h, while 
attempting to leap to an adjacent leaf, the pair fell into a 
spider web. They remained in amplexus as the female tried 
to break free from the web, unsuccessfully. After 60 min 
observing the individuals in the spider web, we collected 
and placed them in a transparent plastic bag with damp 
leaves from the surrounding vegetation. On the next day, 
we found an egg clutch deposited inside the collection bag, 
attached to the upper side of a leaf, with the female not 
brooding the eggs. Finally, in the fourth event on 3 Decem- 
ber 2018, we collected an amplectant pair that we placed in 
a transparent plastic bag for direct observations. The col- 
lected female showed the same behaviours as the females 
observed in nature, with the amplexus lasting ~ 100 min. 
Clutches observed after oviposition and care had a mean of 
14 eggs (13-16), measured 1.57 cm in diameter (1.30-2.09), 
and 4.5 mm in thickness (3-6). Female egg brooding lasted 
for a mean of 110 mins (range = 99-120 mins, N = 3). 

Although we observed some males calling during court- 
ship and amplexus, we were not able to record them. To the 
ear, the calls emitted during these events appeared identi- 
cal to the advertisement call. We checked the clutches of 
the aforementioned events for two days following oviposi- 
tion. We did not observe individuals close to them, neither 
at day nor at night. On the second night, we found that all 
eggs from all these clutches were missing, purportedly pre- 
dated. We did not observe aggressive behaviours between 
conspecific males or predation events involving individu- 
als or clutches. 

In the analysis of abiotic variables affecting reproductive 
activity, we selected the additive model without interaction 
for number of calling males and simple models for number 
of sampled individuals. The number of calling males was 
positively correlated with humidity (X? = 10.6828, P = 
0.0010813, N = 13) and precipitation (X? = 13.8005, P = 


Reproduction and phylogenetic relationships of Vitreorana baliomma 


<0.001, N = 13) and negatively correlated with temperature 
(X? = 6.8499, P = 0.0088646, N = 13) (Fig. 3). The number 
of sampled individuals was positively correlated with hu- 
midity (X? = 8.2314, P = 0.004117, N = 13) (Fig. 4) and not 
related to precipitation (X? = 1.4078, P = 0.2354, N = 13) or 
temperature (X? = 3.1656, P = 0.07521, N = 13). 


Microhabitats of adults and oviposition sites 


We found all adults on the upper faces of green smooth 
leaves (N = 24). The sizes of used leaves ranged from 6.6 x 


2.2 to 39.9 X 8.5 cm. Most adults used the left side of streams 
(N = 23), and were perched on leaves positioned above the 
streambed (N = 14) (Table 1). We found all clutches depos- 
ited on green smooth leaves (N = 19), ranging in size from 
6.0 X 2.4 to 61.3 X 28.6 cm, more frequently on the upper 
face (N= 15) and on the left side of the stream (N = 14), 
above the streambed (N = 15) (Table 1). 

The microhabitats used by adults and oviposition sites 
were significantly different in vertical height (t = 4.2009, 
P = <0.001, N = 43) and horizontal distance (W = 300, P = 
0.02703, N = 43), but not in leaf size (W = 188, P = 0.334, 
N = 43), canopy openness (W = 96, P = 0.3602, N = 32) 


Figure 2. Reproductive behaviour in Vitreorana baliomma. While in amplexus, pairs move over leaves and branches (A), change the 
leaf face (B), and perform 360 and 180° rotations on leaves (C). After oviposition, the female will brood her egg clutch, elevating the 


posterior part of her body (D). Photographs depict various pairs. 
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Table 1. Variables of microhabitats used by adults and oviposition sites of Vitreorana baliomma at two streams in the Atlantic Forest 
in the south of the state of Bahia, Brazil. * 20 adults measured, ** 12 clutches measured. 


Stream side 


Horizontal Leaf surface Position at stream Canopy open- 
Height f Leaf area 5: 
fem) distance Upper Lower lea?) Left Right Center Streambed Bank Outside ness 
iji (em) 06) (6) O (%) A A A (%) bank (%) (6) 
Adults 142+43.5 12.3421 100 0 121+76.2 95 5 0 62 33 5 10.5+3.11 
(N=24) (60.7-231.4) (0-72) (11.4-266.2) (3.72-15.24)* 
Clutches 79452 3.2+12.9 79 21 236+350.5 74 21 5 79 21 0 13.75+8.13 
(N=19) (6-219) (0-58) (11.3-1376.2) (6.58-37.89)** 
Locality 
e BA001 
30 30 30 4 Mascote 
20 20 a 


# calling males 


26 27 
Average temperature (°C) 


80 84 28 29 


Average humidity (%) 


0 10 20 30 88 92 25 


Daily rainfall (mm) 


40 50 


Figure 3. Correlation between the number of calling males of Vitreorana baliomma and daily rainfall (X? = 13.8005, P < 0.001), humidity 
(X? = 10.6828, P = 0.0010813), and temperature (x? = 6.8499, P = 0.0088646). The grey area represents the 95-96 confidence interval. 


(Fig. 5), position at stream (X? = 2.4087, P = 0.2999, N = 
43), stream side (X? = 4.4682, P = 0.1071, N = 43), and leaf 


surface (X? = 3.355, P = 0.067, N = 43) (Fig. 6). 
10 


Clutch characterization 


Mode clutch size was 14 + 4 eggs (7-22). The mean dia- 
meter and thickness of egg capsules were 2.49 + 0.75 cm? 
(1.31-4.32) and 6.7 + 2.8 mm (2-13), respectively (N = 18). 
Clutches had embryos between stages >10 and 25 accord- 
ing to GOSNER (1960), where >10 was the stage immediate- 
ly after oviposition and 25 the stage of hatching embryos. 
We found two differently aged clutches on the same leaf 
(stages 19 and 21). All clutches exhibited a jelly-rich matrix 
surrounding the embryos. At stages 10-16, we observed 
yellowish cream coloured embryos (slightly greenish im- 


# sampled adults 


5 pene i mediately after deposition) (Fig. 7A). At stages 19-22, we 
4 Mascote observed whitish or slightly greyish embryos, with a yel- 
lowish cream vitellus (Fig. 7B). Late-stage embryos (23-25) 

80 84 88 92 


Average humidity (%) 
Figure 4. Correlation between number of sampled adults of 


Vitreorana baliomma and humidity (x? = 8.2314, P = 0.004117). 
The grey area represents the 95-96 confidence interval. 
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were light brown to grey with reddish hues and a reduced 
yellowish cream vitellus (Fig. 7C). In early-stage clutches, 
eggs were in a monolayer array (Fig. 7A). However, this 
changed over time (Figs 7B, C), with most clutches increas- 
ing in thickness during development. 


We selected simple models for the analysis of both 
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clutch structure and microhabitat affecting clutch thick- 


ness. Clutch thickness was positively correlated with devel- 
opmental stage (F = 8.6555, P = 0.009568, N = 18) (Fig. 8). 
However, clutch thickness was not related to the number of 
eggs in the clutch (F = 2.2403, P = 0.1539, N = 18) or clutch 
diameter (F = 0.756, P = 0.3974, N = 18). As for the effects of 
microhabitat data, clutch thickness was positively correlat- 
ed with leaf size (x7 = 5.1464, P = 0.02329, N = 18) (Fig. 8), 
whereas vertical height (x7 = 0.2118, P = 0.6453, N = 18) and 
leaf surface (x? = 1.4638, P = 0.2263, N = 18) did not influ- 


ence clutch thickness. 


Height (cm) 
Š 


100 


Phylogenetic analyses 


Our final dataset included 3176 characters (12S: 984 aligned 
positions, 18% of terminals; 16S: 883, 21%; COI: 705, 89%; 
POMC: 604, 69%). ModelFinder selected the follow- 
ing data partitions and models of nucleotide evolution: 
TIM2+FU+I+G4 for 12S and 16S, TIM2e+I+G4 for COI 
first codon and POMC first codon, TN+FU+I for COI 
second codon and POMC second and third codons, and 
TN+FU+G4 for COI third codon. ML searches found a 
single most likely tree (Ln = -18153.711; Fig. 9), which re- 


covered a monophyletic Vitreorana (BS = 80) within 
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Figure 5. Comparison between quantitative variables of the microhabitat used by adults and oviposition sites. Median (horizontal 
line), quartiles (boxes), extremes (vertical lines), and outliers (points) are shown. Significant values are pointed out as * P < 0.05 and 
e P< 0.001. 
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Figure 6. Comparison between qualitative variables of the microhabitat used by adults and oviposition sites. 
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Centroleninae (BS = 70). All species of Vitreorana with 
more than one sampled specimen are monophyletic ex- 
cept for V. uranoscopa, which is paraphyletic with regard to 
V. franciscana. Within Vitreorana, the clade formed by the 
Guiana Shield/Amazonian species V. helenae and V. ritae 
(BS = 100) is sister to all other Vitreorana species (BS < 50). 
The two Cordillera de la Costa species (V. antisthenesi and 
V. castroviejoi) form a clade (BS = 100) sister to all the other 
species. The Atlantic Forest and Cerrado species of Vitreo- 
rana are not monophyletic because V. baliomma is sister to 
V. gorzulae from the Guiana Shield (BS < 50). 


Discussion 
Reproductive behaviour and parental care 


The briefness of female brooding after amplexus is prob- 
ably key to understanding why this behaviour went un- 
noticed in the family Centrolenidae until now. Through 
detailed field studies and careful observations during ovi- 
position, our study was the first to detected parental care 
in an Atlantic Forest glassfrog, broadening the knowledge 
about the reproductive biology for the group. Vitreorana 
baliomma females perform brief brooding just after ovipo- 
sition. Although we have not quantified the effects of this 
behaviour on embryo survival, it is likely to have a similar 


impact to what DELIA et al. (2017) found for Cochranella 
granulosa and Teratohyla pulverata, two species of close- 
ly related genera (GUAYASAMIN et al. 2009), where brief 
brooding functions to hydrate eggs, protecting embryos 
from dehydration and early predation. Thereafter, clutches 
will depend on water from rain and ambient humidity to 
maintain hydration levels and clutch thickness (DELIA et 
al. 2017, 2020). This discovery in V. baliomma reinforces the 
pattern of parental care evolution hypothesized for centro- 
lenids (DELIA et al. 2017) and suggests that the behaviour 
may be widespread in Vitreorana and Cochranellini spe- 
cies. Nevertheless, more studies and experiments are nec- 
essary to determine how brooding impact the clutches. 
During brooding, we observed females elevating the 
posterior part of their body. This movement may serve to 
trap humidity (considering that the belly of anurans has 
a greater absorption capacity), reposition the body above 
the clutch, or aid egg fertilization. Other centrolenid spe- 
cies with brief maternal care brooded for similar periods 
of time. For example, the female of V. ritae brooded for 80 
min (DELIA et al. 2017). However, variation is large, from 10 
to 180 min, in species of the genera Centrolene, Espadara- 
na, Nymphargus, and Teratohyla (JACOBSON 1985, GUE- 
VARA-MOLINA & VARGAS-SALINAS 2014, D1Az-RICAURTE 
et al. 2016, ARCILA-PEREZ et al. 2017, BASTO-RIASCOS et 
al. 2017, R1os-SoTo et al. 2017). Some species rotated their 


ay 


Figure 7. Egg clutches of Vitreorana baliomma at early (A), intermediate (B), and late (C) stages of development. 
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eggs while caring, but ours is the first record of this body 
elevation behaviour. 

Our discovery of female brooding behaviour in V. bali- 
omma contrasts with the only two studies (GouvEIA et al. 
2012, LISBOA et al. 2019) that provided data on reproduc- 
tive aspects of V. baliomma, but did not document paren- 
tal care in this species. LISBOA et al. (2019) did not observe 
oviposition events, which significantly reduces the chances 
of observing female brooding behaviour. Gouveta et al. 
(2012), based on a single amplectant pair, described that 
the female stayed near the clutch until dawn, a behaviour 
we did not observe, without mentioning interactions be- 
tween the female and her egg clutch. This could indicate 
that V. baliomma regulates its maternal care based on envi- 
ronmental factors such as ambient moisture, daily rainfall, 
or both. This would be similar to the variation in pater- 
nal care observed in Hyalinobatrachium (VOCKENHUBER et 
al. 2008; DELIA et al. 2013) or attributed to some degree 
of behavioural plasticity (WELLS 1981). Another possibility 
is that maternal brooding is subject to variation between 
certain populations, which means that this character could 
have a phylogeographic signal function, possibly with a 
genetic component influencing this behaviour (BUBAC et 
al. 2020). Finally, we cannot rule out that these authors 
disturbed the females when taking measurements of the 
eggs immediately after oviposition as stated in GOUVEIA 
et al. (2012: 18). Individuals engaged in caring are normal- 


ly more risk-tolerant, as females staying near their clutch- 
es will surely potentially attract the attention of predators 
(BRAVO-VALENCIA & DELIA 2016, VALENCIA-AGUILAR et 
al. 2020b). It is important to test if V. baliomma exhibits 
variation in parental care and, if so, which variables could 
be associated with it. 

The general sequence of observed mating behaviours 
is shared with all other species of glassfrogs, with varia- 
tion being associated with details regarding aspects such 
as call and oviposition sites, caring sex, and extension 
and type of parental care (GREER & WELLS 1980, JACOB- 
SON 1985, VOCKENHUBER et al. 2008, VARGAS-SALINAS et 
al. 2014, ARCILA-PEREZ et al. 2017, BASTO-RIASCos et al. 
2017, R1os-SoTo et al. 2017, Diaz-RICAURTE et al. 2019). 
The behaviours observed during the amplexus suggest 
that the female selects the best place for depositing the egg 
clutch. While in amplexus, the female does not only move 
through the foliage of the plant in which the amplexus 
started, but also through that of other plants nearby and 
frequently switches between leaf faces. This face-changing 
behaviour explains why we found 21% of clutches on the 
undersides of leaves. This could be random variation, but 
we argue that it could also result from a trade-off, where 
eggs on the upper face will be exposed to rain and so can 
benefit from a non-parental source of hydration, while 
those on the lower face may be better protected from em- 
bryo loss. 
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Figure 8. Correlation between clutch thickness of Vitreorana baliomma and developmental stage (F = 8.6555, P = 0.009568) and leaf 
size (X? = 5.1464, P = 0.02329). The grey area represents the 95-96 confidence interval. 
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Amplectant pairs moved about for more than an hour 
before oviposition. This is compatible with the observa- 
tion that species with maternal care move across greater 
distances, investing more time in amplexus than those 
with paternal care (JACOBSON 1985). This might be related 
to territoriality, inasmuch as in species with paternal care, 
males might choose a place that is good for both signal- 
ling for mates and embryonic development. As these sites 
are probably available only to a limited extent, males are 
expected to display more aggressive territorial behaviours 
(McDIARMID & ADLER 1974, RESTREPO & NARANJO 1999, 
CARDOZO-URDANETA & SENARIS 2012, but see NORONHA 
& RODRIGUES 2018) and do not move far away from their 
territory. On the other hand, in species with maternal care, 
males might choose the best site to enhance their acoustic 
(or other) signals, while females must choose a good place 
for embryo development (ToUCHON 2012). 

As expected, we found a positive relationship be- 
tween average relative humidity and daily rainfall with the 
number of calling males. This relationship points to the 
importance of precipitation for breeding activity in V. bali- 
omma. Notwithstanding, we found a negative relationship 
between average temperature and the number of calling 
males. This probably results from an inverse correlation of 
temperature and daily rainfall. Sampled adults presented 
a relationship with relative humidity too, with more in- 
dividuals found at times of higher humidity. Despite our 
sample period having a relatively low frequency and a non- 
standardized duration of sampling, our results were com- 
patible with other studies analyzing the reproductive ac- 
tivity in centrolenids. Precipitation is a good predictor of 
calling activity in the Centrolenidae (VoCKENHUBER et al. 
2008, Basto-Riascos et al. 2017, Rros-SoTo et al. 2017), 
and Vitreorana species call, almost exclusively, during and 
after rainy days (CARDOSO & HADDAD 1992, CANELAS & 
BERTOLUCI 2007). Also, activity of male and female centro- 
lenids tends to be higher in conditions of more rainfall and 
lower temperature (BAsSTO-R1Ascos et al. 2017, Rios-SoTo 
et al. 2017). Vocalizing is energetically demanding for an- 
urans, causing a significant increase in their metabolic rate 
(GERHARDT 1994), so that it is expected that species call 
just when the environmental conditions are good for do- 
ing so. In fact, during the unusually long dry periods (i.e., 
more than three consecutive weeks of almost zero precipi- 
tation) during our sampling, the entire population became 
inactive, and no males were heard, or individuals observed. 


Microhabitats of adults and oviposition sites 


Our results on variables associated with the oviposition site 
indicate a preference for the upper sides of large smooth 
green leaves, relatively close to the substrate, above water, 
and on the left side of streams. Meanwhile, adults also ap- 
pear to prefer the upper side of smooth green leaves, above 
water, and on the left side of streams, but in contrast to ovi- 
position sites exhibit a preference for greater heights above 
the ground, larger horizontal distances from water and, ap- 


parently, smaller leaves. Analyses comparing the micro- 
habitats of adults and clutch deposition sites revealed sig- 
nificant differences only for vertical height and horizontal 
distance. However, our sample size is too limited for egg 
clutches to comprehensively confirm patterns of micro- 
habitat use, especially for the preference of stream side, 
leaf size, and canopy openness. Oviposition sites located 
at lower heights can result in higher survival probabilities 
(JACOBSON 1985), perhaps due to higher humidity levels 
closer to the ground (CaARDOso & HADDAD 1992), which 
could decrease dehydration. Furthermore, these sites could 
benefit tadpole hatching by decreasing the drop distance 
to the water and the influence of wind. On the other hand, 
because most adults that we observed were calling males, 
we hypothesize that their greater perch height, compared 
to clutches, enhances the dissemination of acoustic signal- 
ling (GREER & WELLS 1980, but see VALENCIA-AGUILAR et 
al. 2012, VARGAS-SALINAS et al. 2014). 

We found all clutches either directly above the stream- 
bed or, less frequently, on the stream bank (in the area that 
will be flooded after strong rainfalls). Considering that 
tadpoles need to fall into the water after hatching to con- 
tinue their development, oviposition sites above streams 
should be strongly favoured by females (CARDOZO-UR- 
DANETA & SENARIS 2012). Meanwhile, we argue that the 
selection of males’ calling sites distance from the water is 
less stringent, allowing males to vocalize from places out- 
side the flooding area, which are more densely covered 
by vegetation and so have a reduced predation risk. Male 
choice for greater heights could be a response too for plac- 
es more covered, once dense vegetation reduces propaga- 
tion of calls (WELLS & SCHWARTZ 1982, VOCKENHUBER et 
al. 2008). 

Although adult microhabitat preferences and oviposi- 
tion site requirements do not significantly differ as to leaf 
size used, clutches were more commonly found on larg- 
er leaves. The use of large leaves might be associated with 
male territoriality and prolonged care in some Hyalino- 
batrachium species, because males can then attend many 
clutches at the same time (VOCKENHUBER et al. 2008, VA- 
LENCIA-AGUILAR et al. 2012). We have not observed male 
territoriality and prolonged care in V. baliomma, but our 
linear model shows a correlation of leaf size with clutch 
thickness. Globular clutches (like in species of Cochranella, 
Espadarana, Nymphargus, and Vitreorana) can retain larg- 
er amounts of water, requiring an initial bout of maternal 
brooding and subsequent rainfall to maintain hydration 
levels (DELIA et al. 2020). We hypothesize that larger leaves 
may accumulate more water on their surfaces, aiding in 
clutch hydration. Thus, water on leaves after rainfall might 
help to ensure embryo survival. 

Clutches and adults did not significantly differ in oc- 
cupying leaf face, although adults were found just on the 
upper side and clutches on both faces. Species with brief 
parental care show tendencies to deposit clutches on the 
upper faces of leaves or other exposed locations (DELIA et 
al. 2020). Although oviposition strategies may be flexible 
with females adjusting their strategies according to ecolo- 
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gical conditions (MARSH & BORRELL 2001), we found that 
V. baliomma females appear to have a preference for de- 
positing their clutches on the upper side, but do not do so 
exclusively. As previously discussed, oviposition on both 
the upper and lower sides of leaves may be related to a 
trade-off between hydration from rainwater on the upper 
side and protection from predators and the washing-away 
of eggs by heavy rainfall and/or shielding from UV light on 
the lower side (BLAUSTEIN & BELDEN 2003, VONESH 2005). 
Leaf-side use for oviposition is considered a taxonomical- 
ly informative trait in the Centrolenidae (e.g., Ruiz-Car- 
RANZA & LYNCH 1991, GUAYASAMIN et al. 2009). Howev- 
er, field observations indicate that this character is more 
intraspecifically variable than previously thought (DELIA 
et al. 2010, BRAVO-VALENCIA & DELIA 2016, ARCILA-PE- 
REZ et al. 2017, MAYNARD et al. 2020). Taking this variation 
into account may have an important impact on our under- 
standing of the evolution of this strategy. 

We suspect that a smooth leaf surface facilitates that 
tadpoles more readily slide and fall into the stream below, 
while the apparent preference of the left side of the stream 
could be related to slope level, which, by pure chance, was 
similar in the two sampled streams, causing the left side to 
receive more water. We did not observe a major difference 
in the vegetation structure between the both stream sides. 
It is likely that, as females do not afford their clutches pro- 
longed care, the selective pressures might act in a way that 
leads to selecting oviposition sites at which clutches will 
be better protected from dehydration. Yet, oviposition site 
choice is a complex process, where females have to con- 
sider the maximization of embryo survival, maximization 
of maternal survival, optimization of offspring phenotype, 
proximity to suitable habitat for offspring, maintenance 
of natal philopatry, or, in an indirect manner, mate choice 
(REFSNIDER & JANZEN 2010). It is therefore crucial to test 
the hypotheses for decisions that rule oviposition site se- 
lection, the impact of female choice (i.e., upper vs. lower 
side of leaves, large vs. small leaves) on offspring survival 
and include observations of other populations and other 
Vitreorana species to improve our understanding of micro- 
habitat use in this group. 


Egg clutches 


Egg and embryo colours, number of eggs, and early stages 
of clutches of V. baliomma are similar to the descriptions of 
GouvEIA et al. (2012), although we noted temporal modi- 
fications to the egg clutches during their development. In 
their early stages, clutches are in a monolayer array; how- 
ever, throughout embryogenesis, they take on a more com- 
plex tridimensional structure due to increased egg and 
egg-jelly hydration. We also found a relationship between 
developmental stage and clutch thickness. As is the case 
with clutch structure, this results from the addition of en- 
vironmental water, as jelly-rich clutches are fast to absorb 
water, store much water, and take more time to dehydrate 
(DELIA et al. 2013, 2020). 
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Phylogenetic relationships 


We here present the first phylogenetic analysis including 
data of Vitreorana baliomma. Maximum likelihood analy- 
sis corroborates the monophyly of Vitreorana, with mod- 
erate bootstrap support. The phylogenetic relationships 
within Vitreorana are different in several important aspects 
from previous well-sampled phylogenies of the genus. Our 
results indicate that the position of V. franciscana within 
the genus is particularly problematic. This species is just 
represented by DNA sequences of POMC (Santana et al. 
2015), a nuclear protein coding gene without sufficient var- 
iation to resolve the relationships between the species of 
Vitreorana. This sampling bias most likely explains the in- 
ferred phylogenetic position of our two samples of V. fran- 
ciscana, which are nested within specimens of V. urano- 
scopa on a conspicuously long branch. Thus, an important 
next step to improve our knowledge of the evolutionary 
relationships of Vitreorana requires gathering more phy- 
logenetically informative DNA sequences and genes, es- 
pecially for V. franciscana, such as mitochondrial genes 
and faster-evolving nuclear markers. The sister relation- 
ship of V. baliomma and V. gorzulae is remarkable. From a 
morphological perspective, both species share completely 
transparent parietal and urinary bladder peritonea, white 
peritonea covering the heart, bulbous liver, stomach, and 
intestines, green bones in life, and dorsum lavender in pre- 
servative (CASTROVIEJO-FISHER et al. 2008, 2009, GUA- 
YASAMIN et al. 2009, PONTES et al. 2014). Considering that 
all other species of Atlantic Forest and Cerrado Vitreo- 
rana (i.e., V. eurygnatha, V. franciscana, V. parvula, and 
V. uranoscopa) share the presence of white pigment on the 
urinary bladder peritoneum, which is absent from all oth- 
er known glassfrogs, the recovered phylogenetic position 
of V. baliomma supports this character as an unambigu- 
ous and non-reversed synapomorphy of all other Atlantic 
Forest species. From a biogeographic perspective, the sister 
relationship of V. baliomma and V. gorzulae indicates a sce- 
nario of multiple dispersal/vicariant events between Ama- 
zonia and the Atlantic Forest, which deviates from the sin- 
gle dispersal into the Atlantic Forest from Amazonia pro- 
posed by CasTROVIEJO-FISHER et al. (2014). We highlight 
that this new scenario should be tested with the addition of 
more characters and other potentially closely related spe- 
cies from the Guiana Shield such as “Cochranella” duidaea- 
na and “Cochranella” riveroi. 
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